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SUMMARY 


The relationship between Joule heating, diffusion fluxes, and friction forces has been studied for both 
total and electron thermal energy equations, using general expressions for multicomponent diffusion in 
two-temperature plasmas with the velocity dependent Lorentz force acting on charged species in a 
magnetic field. It is shown that the derivation of Joule heating terms requires both diffusion fluxes and 
friction between species which represents the resistance experienced by the species moving at different 
relative velocities. It is also shown that the familiar Joule heating term in the electron thermal energy 
equation includes artificial effects produced by switching the convective velocity from the species 
velocity to the mass-weighted velocity, and thus should not be ignored even when there is no net 
energy dissipation. 


1. INTRODUCTION 


Electrical current in a plasma that is a mixture of different types of ions, electrons, and neutral species 
produces energy dissipation, i.e.. Joule heating, which directly contributes to the thermal energy or 
temperature of plasmas. This term is usually represented as E in energy equations, where and E 
represent electrical current density and electric field, respectively. In the presence of a magnetic field 
B, the Joule heating term includes the effect of the Lorentz force with results J^ (E + u x B), where u 
is the mass-averaged velocity of the plasma. 

Due to the small mass of electrons, collisional energy transfer is inefficient between free electrons and 
heavy particles. As a result, persistent temperature differences often exist between free electrons and 
heavy particles, and this thermal nonequilibrium requires a separate energy equation for free electrons 
to be included in the analysis. Joule heating in the electron energy equation presented in textbooks 
(refs. 1, 2) takes the form of + u x B), where q r J, =p, (u, - u), and u, are respectively the 

electrical charge per unit mass, diffusive mass flux, and average velocity of species /, and subscript e 
denotes free electrons. Since qj t is the electrical charge carried by diffusion flux, q ( J,-(E + u x B) 
appears to represent the work done by the motion of electrical charge through E and B. 



A question arises, however, regarding the origin of these terms. Forces acting on the charged species 
due to E and B are body forces, and terms representing the work done by body forces should not 
appear in thermal energy equations. This is precisely the reason that the work done by gravitational 
forces does not appear in thermal energy equations. It simply drops out when the mechanical energy 
equation for potential and kinetic energies, i.e., Bernoulli’s equation, is subtracted from the total 
energy equation to obtain the thermal energy equation. Joule heating, however, undeniably contributes 
to the thermal energy, and its appearance in the thermal energy equation is therefore physically correct, 
contradicting the fact that the work done by body forces such as electromagnetic forces should not 
appear. 

In the electron thermal energy equation widely used in high-temperature gas dynamics (ref. 3), the 
Joule heating term is simply treated as work done by body force (electric field) and added as a source 
term in the electron thermal energy equation. Also in a popular magnetohydrodynamic (MHD) 
description (ref. 1), the ^/(E + u x B) term is brought into the electron total energy equation as the 
work done by body forces, and kept in the thermal energy equation by simply ignoring the electron 
inertia rather than subtracting the mechanical energy from the total energy. Obviously, both of these 
approaches are unsatisfactory for the reasons explained above. 

In contrast, multifluid descriptions (refs. 4, 5) treat the energy dissipation caused by the velocity 
difference between electrons and ions as the sole source of Joule heating instead of the work done by 
body forces. The physical mechanism of Joule heating is correctly represented in those descriptions, 
since Joule heating is not related to the work done by body forces, but to the frictional dissipation, as 
shown in this paper. In order to represent Joule heating as J^, E and B, however, it is necessary to 
include additional terms involving diffusion fluxes brought into the system by switching convective 
velocity of the species from u, to u. It is not, in general, possible to simply replace the frictional 
dissipation term to J ? -(E + u x B). 

Another question arises regarding the qj r i E + u x B) term in the electron thermal energy equation, 
in particular, for the special case of zero electrical current density, i.e., = I, <?, J, = 0, which is 

commonly referred to as “ambipolar diffusion.” During the recent development of formulations for 
multicomponent diffusion in plasmas (refs. 6—8), it has been shown that concentration gradients of the 
charged species in plasmas result in nonzero E, even when J v = 0. Here, the Joule heating term in the 
total thermal energy equation is, of course, zero, since J v = 0. The ^/(E + u x B) term in the 
electron thermal energy equation, however, appears to be nonzero even in this special case. This 
counterintuitive behavior can also be explained by examining the relationship between diffusion and 
energy dissipation due to friction between species. 

This paper presents a study of the relationships between Joule heating, diffusion fluxes, and friction 
between species for total and electron thermal energy equations. It is shown that Joule heating results 
from the combination of both diffusion fluxes and friction between species, which represents the 
resistance experienced by the species moving at different relative velocities. Use of the hydrodynamic 
theory of diffusion (refs. 8—10), instead of the standard kinetic theory (refs. 11-16), has provided a 
clear and simple physical interpretation of complex phenomena of multicomponent diffusion and 
frictional interaction between species in two-temperature plasmas with nonzero electrical current and 
electromagnetic fields. Section 2 presents a derivation of the total and electron thermal energy 
equations from the species thermal energy equation. The Joule heating terms in total and electron 
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thermal energy equations are constructed from friction forces and diffusion fluxes brought into the 
system by artificially switching convective velocities from u, to u. In addition to the real energy 
dissipation (Joule heating), the q t J f -(E + u x B) term in the electron thermal energy equation therefore 
includes this artificial component, and should not be ignored even when = 0. Section 3 presents 
physical interpretations regarding the relationships between Joule heating, diffusion fluxes, and 
friction between species. Section 4 contains a few concluding remarks. 


2. JOULE HEATING IN ENERGY EQUATIONS 


Chemical reactions and viscous dissipation are not included for simplicity in the present analysis, since 
they are not directly related to Joule heating and diffusion fluxes. The thermal energy equation for each 
species i is then given by (ref. 13) 


d(P, g ,) 

dt 


+ V(p,q u,) = 


- V • q, - pV u, + £ Q n + X £„(“ 

/ i 


U,) F„ 


( 1 ) 


where p ( and p, are respectively the partial mass density and partial pressure of species /; e i (T l ) is the 
partial specific' internal energy of species i at temperature T;, and q, represents the conduction heat flux 
vector Q = -Q is the energy exchange between species i and j, and F, ; = - F ;( is the mean “friction 
force per unit volume of species j on species i with they summation extending over N components in 
the mixture. Note that Q„ = 0 and F„ = 0. General expressions for heat flux vector q, have apparently 
not been derived in full generality, but results are available in particular cases (refs. 2, 4, 14-16). 
Nevertheless, it is not necessary to know explicit forms of q, and Q tJ for current purposes, (u^ - u,> „ 
represents the energy dissipation by the friction between species i and j, and is its fraction into 
species i which is given by (ref. 13) 


m, + m i 

where m, is the mass of a single particle of species /. F, y is of the general form (ref. 8) 

F,y = °v( u / " u,)+y„ (: 

Y,i = Pa • V In 7’ / - ^ • V In T, (* 

where a and p are frictional and thermophoretic force coefficients respectively, which now become 
tensors due to the presence of B (ref. 4). a tJ is related to binary diffusion coefficient tensors D, y as 

(ref. 8) 

D. = P=rf*;' ( 
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where z, = p j Ip. Derivation of equation (1) involves all body forces, including electromagnetic forces 
(ref. 13). 

Although energy dissipation due to friction forces between species is represented by (u ; - u,) F y , the 
Joule heating term represented by J^, E, and B does not explicitly show up in equation ( 1 ) as 
expected, for Joule heating is not related to the work done by the body force. Joule heating repre- 
sented as J V (E + u x B) or + u x B) appears only in single-fluid or MHD descriptions as a 

result of switching the convective velocity from u, to u, as the following development shows. First, 
equation (1) is rewritten with u as the convective velocity, with results 

+ V • (p,e,u ) = -V • (q, + h,J,)~ P y ■ u + — ■ Vp, + £ Q, t .+ XC/( U F „ (6) 

ot P, , j 

where h t = e, + p, Ip , is the partial specific enthalpy of species 


A. Total Thermal Energy Equation 

The total thermal energy equation is obtained by summing equation (6) over all species, with the result 

^^ + V-(peu) = -V •J(q, + A,J,)-pV>u+J- -Vp, + J,C,( U ( 7 ) 

ot i , P, // 

where e = (l/p)I,p,e, is the total thermal energy, p = S, p, is the total pressure, and use has been 
made of 1 , 0, = 0. Derivation of the Joule heating term consisting of J v , E, and B requires both 
I, (1/p,) J ,-Vp, and Z, y ^ (u y - u,)-F , 7 in equation (7), as shown below. 

In the limit of large friction between species, in which diffusional behavior results, the partial pressure 
gradient is given by (ref. 8) 


Yp. 

p, 


Du _ 1 

+ F, + — 

Dt p, 



(B) 


where D/Dt = d/dt + u V, and F, is the body force per unit mass acting on species i which is taken to 
be of the form 


F = g + # + u,xB) 


(9) 


where g is the acceleration of gravity. 

Using equation (8), it is shown that 

I- Vp, = Xv.J. (E+u ; xB)+£(u,-n) F„ (10) 

l Pi I </ 
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Use has been made of X, J, = 0. It is also shown that 


J, (u,xB) = J, 


— x B + 
vA 


u x B j 


= J, (uxB) 


Frictional terms are combined together with the result 


(ID 


X C„(u , - u , ) • F, + 5> , - u ) • F ( . = X 

// »/ 1 / 


'", u , + w , u , 


if 


Observe that X jy = 0, since A tj = - A p . From equations ( 10)-( 12), it is shown that 


I- ■ ?/>, + IC,<U u,) F = J, - (E+ u x B) 

/ Pi ij 


( 12 ) 


(13) 


The total thermal energy equation then becomes 

+ v • (pcu) = -V • ^(q, + /i,J,)- pV • u + J v - (E+ u x B) (14) 

dt , 

Equation (14) is the standard total thermal energy equation which can be found in textbooks 
(refs. 1, 2). As shown above, construction of Joule heating requires both (1) frictional forces that 
represent the collisional coupling between dissimilar species moving at different relative velocities 
and (2) diffusion fluxes brought into the equation by switching the convective velocity from species 
velocity u, to the mass-averaged velocity u. 


B. Electron Thermal Energy Equation 


The electron thermal energy equation is obtained by simply replacing the subscript i with e in 
equation (6) with the result 


+ V • (p e,u ) = -V • (q 4 . + h, J , ) - p,V • u + ^ • Vp t . + £ Q,, + X C,( u , " u ) ' ' F « 
dt Pc 


(15) 


Note that equation (15) contains (1 /p f ) J and X ( C, a (u, - u^-F^, instead of the more familiar 
q J -(E + u x B) term as in equation (5.2) of section IV in reference 2. As in the case of the total 
thermal energy equation, the ^J f (E + u x B) term is obtained by combining ( 1/p,,) J and 
£. (u, - u f )-F w and by further simplifications as shown below. 


Replacing Vp,, in equation (15) using equation (8) would not make equation (15) simpler, since more 
terms would be brought in. Nevertheless, it is possible to seek further simplification by exploiting the 
fact that the electron mass m e is very small compared to the masses of heavy particles. This can be 
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done by neglecting terms of order £ 2 = m f in comparison to terms of order unity. Equation (8) for 
electrons then becomes 


^ = 9 f (E+u f xB)+- JX (16) 

P, P, ■ 


Using equation ( 1 1), it is shown that 

— • Vp, = q e J , • ( E + u x B ) + • X P,, ( 1 7 ) 

Pe Pe ’ 

F ( ., is not ignored, since a n is of order £ = ^rn L , (refs. 7, 8). Using equation (17) and £ * 1, 
equation (15) then becomes 

^^ + V (p A u)=-V (q , + y < )- ft V-uta-(E + «*B) 

* ( 18 ) 
+ + 

/ / 

Equation (18) becomes essentially the same as equation (5.2) of section IV in reference 2, when 
Z, (I/#) J,-*F ri is ignored, which can be done by neglecting either J, or F w , although there is no 
apparent justification for doing so. 

As briefly described in section 1, the ^-(E + u x B) term, originally brought into the system as the 
work done by body forces, is kept in the electron thermal energy equation with u e as the convective 
velocity (eq. (5. 196) in ref. 1). This is done by ignoring the inertia of electrons instead of following 
the more general practice of subtracting the mechanical energy from the total energy, which would, of 
course, eliminate the qj/i E + u x B) term in the thermal energy equation. However, equation (17) 
illustrates that ignoring electron inertia also means that q t J/(E + u x B) drops out anyway from 
equation (5. 196) of reference 1 by canceling with the (u ( . - tO-Vp, term, which is also present in 
equation (5.196) of reference 1. (All frictional interactions (F„) are also neglected when electron 
inertia is neglected in reference 1 , which would result in an ideal MHD description without net energy 
dissipation as shown below.) By using u as the convective velocity, the q t J f -(E + u x B) term is, of 
course, brought back into the electron thermal energy equation, as the preceding development shows. 

Further simplification could be explored by neglecting terms of order £ = as well. Equation (8) 
then becomes 


— = E+u,xB+ — Yr,,, (19) 

PA- ‘ PA , 

Observe that y a = /3 ie -VlnT tJ , since (i u . is order unity, and j 9 ei is order e ■ (ref. 17). Using 
equation (19), equation (15) is rewritten, with results 
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i^el + v . {p e u) = -V . (q + h e J e ) - p e v ■ u + q e s e • (E + u x B) 

dt 


( 20 ) 


+ X Qei + 

f / 

Equation (20) is not much simpler than equation (18), unless thermal diffusion of electrons is 
negligible. Furthermore, this simplification of ignoring terms of order e, would result in an ideal MHD 
description in which resistive effects vanish, since the electrical conductivity of the plasma o is itself 
of order £ 1 (refs. 8, 12). Observe that equation (19) is essentially equivalent to the a limit of 
equation (9.47) of reference 18, which indeed confirms the passage into the realm of ideal MHD. 

(In the present context, however, the term “ideal” does not imply that the plasma flow is isentropic. 
Diffusion and finite-rate chemical reactions are still irreversible processes.) Note that it still is possible 
to keep resistive effects by ignoring only terms of order e 2 , while keeping terms of order £, which is 
precisely the way cf is represented by the binary diffusion coefficients of electrons and ions (ref. 12). 


3. JOULE HEATING IN MHD 


It can easily be explained that the acceleration of the charged species due to electromagnetic forces and 
the resulting increase of the kinetic energy would eventually convert to the increase in random thermal 
velocities through collisions (dissipation) between dissimilar species. Friction between the same 
species is zero (F l( = 0), and applying body forces such as gravitational or electromagnetic field does 
not result in nonzero F„, since all the particles of the same species are subject to the same amount of 
forces and acceleration. Therefore, applying g, E, or B merely changes u ( rather than random thermal 
velocities, and does not contribute to the thermal energy without a proper dissipation mechanism such 
as friction between dissimilar species. In plasmas composed of both positively and negatively charged 
particles, particles with positive charges are accelerated by E and B in the opposite direction ot the 
acceleration of particles with negative charges, increasing velocity differences which in turn produces 
electrical current. This increase in velocity differences is balanced by increased frictional forces repre- 
sented as resistivity. In other words, the “extra” source of acceleration or driving forces due to E and 
B is balanced with “extra” frictional forces, and continued supply of force and energy by E and B 
results in the increase of thermal energy by friction. This clearly shows that the resistance experienced 
by species moving at different velocities is an essential mechanism in Joule heating. 

As shown in the preceding sections, however, diffusion fluxes are also necessary in order to represent 
the Joule heating in terms of E, and B in the energy equations. These diffusion fluxes are brought 
into the system of single-fluid or MHD descriptions, since species in the plasma in these descriptions 
are thought to be moving with the convective velocity of the mass-weighted velocity u of the mixture 
rather than with their own velocities u,. If u is kept as the convective velocity, there is no Joule 
heating term consisting of J„, E, and B in the energy equations as in equation (1), consistent with the 
fact that the work produced by body forces should not appear in thermal energy equations. 

In multifluid descriptions, in which each species has its own momentum and energy equations, Joule 
heating terms do not explicitly show up in the species thermal energy equations. Nevertheless, Joule 
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heating is provided to the system by acceleration terms in the species momentum equation and 
frictional dissipation in the species energy equation, as explained above. The electron energy equation 
in widely used multifluid descriptions (refs. 4, 5) involves the frictional dissipation in the system, 
consistent with the mechanism described above. However, simple conversion of the frictional term to 
the Joule heating as in equation (2. 18) of reference 4 is valid only when Vp e is ignored, as shown in 
equation (17). This can also be easily observed by inspecting equations (2.2e) and (2.18) in 
reference 4. 

In contrast, the increase in differences between species velocities due to the electromagnetic field is not 
represented by the single momentum equation for the whole mixture used in single-fluid or MHD 
descriptions (refs. 1, 2). Diffusion fluxes nevertheless include the increase in species velocity differ- 
ences caused by extra driving forces due to E and B which are, of course, balanced with friction 
between species as described above. The species velocity differences and frictional dissipation eventu- 
ally become the Joule heating term in energy equations, when diffusion fluxes are combined with the 
frictional dissipation terms as shown in the preceding section. 

When J = 0, the Joule heating term in the total thermal energy equation is, of course, zero. However, 
the qjji E + u x B) term in equation (18) or (20) does not vanish even in this special case as the 
question raised in section 1 shows. The fact that concentration gradients of the charged species in 
plasmas result in nonzero E, even when = 0, can easily be illustrated by equation (19). In 
particular, in the absence of magnetic field and temperature gradient, equation (19) simply becomes 


PA 


This has the simple and intuitive interpretation that electrons are in an electrostatic equilibrium with the 
electric field E, which is precisely analogous to hydrostatic equilibrium in a gravitational field. E in 
equation (21) is not induced by charge separation as discussed in a previous study (ref. 19). On the 
contrary, charge neutrality in the plasma is preserved by this E (refs. 6, 7). 

The preceding development shows that q c J ( (E + u x B) is obtained by combining the frictional 
dissipation represented by I, (u, - u,)-F f , with ( 1/p,) J/Vp, which is artificially brought into the 
system by switching the convective velocity from u f to u. Since ( 1/p,) J is in general nonzero, 
the qJ e (E + u x B) term obviously is nonzero as well. In other words, the qJ e -(E + u x B) term 
involves an artificial effect of switching the convective velocity in addition to the energy dissipation by 
friction, and thus should be included in the electron thermal energy equation even when = 0. 

Note that equations (19) and (21) are derived by neglecting the terms of order £ in equation (16). Since 
the frictional force coefficient a ei is also of order £, a a vanishes from the system, which in turn results 
in infinite electron diffusion coefficients (eq. (5)) and electrical conductivity (refs. 8, 12). In other 
words, zero simply deactivates frictional interaction between electrons and other species, and 
frictional dissipation and resistive effects vanish. In this simplification, </,J t ,(E + u x B) merely 
replaces ( 1/p,) J ,-Vp r and no energy is dissipated by this term. 
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When J and m e are not zero, the correct amount of Joule heating is produced by the "extra diffusion 
of charged species induced by “extra” driving forces, as explained above. This can be automatically 
accomplished, when J„ E, and B are consistently determined by solving the generalized Stefan- 
Maxwell equations (ref. 8) involving nonzero m e and Maxwell’s relations in a coupled manner. The 
solution of these equations determines E as well as J,, and this implicitly determines the relation 
between E and J v , i.e., Ohm’s law. This shows that Ohm's law is not an independent constitutive 
relation in multicomponent MHD, and it is not, in general, possible to express J„ directly in terms of 
E, e.g., “generalized” Ohm’s law, without determining J,. 


4. CONCLUDING REMARKS 


A study has been presented of the relationship between diffusion fluxes, friction forces, and Joule 
heating in two-temperature multicomponent MHD, using hydrodynamic theory of diffusion 
(refs. 8-10) which provides a clear physical interpretation. It has been shown that energy produced 
by J flowing through E and B contributes to thermal energy via additional driving forces for diffusion 
by E and B and the resistance experienced between dissimilar species moving at different relative 
velocities rather than the work done by body forces as described in previous studies (refs. 1-3). It is a 
coincidence that the 9( J,-(E + u x B) term brought into the electron thermal energy equation as the 
work done by body forces (refs. 1-3) is identical to the term derived in this paper from frictional 
dissipation and diffusion fluxes. 

Derivations of both total and electron thermal energy equations have also been presented, and it has 
been illustrated that the derivation of Joule heating terms requires both diffusion fluxes and friction 
forces between species. The qJ e -(E + uxB) term in the electron thermal energy equation involves 
the (1/p ) J ,-Vp term, which is artificially brought into the system by switching the convective 
velocity from u, to u. (1/p,) J ,Vp, is not, in general, negligible, and qJ r -( E + u x B) therefore 
should not be ignored, even when = 0. 

Ignoring frictional interaction between species, in particular between electrons and other species, leads 
to equations (19) and (21). Evaluating E using equation (19) or (21) results in an ideal MHD descrip- 
tion and thus should be avoided to keep resistive effects. The proper amount of Joule heating would 
automatically be produced when J„ E, and B are evaluated by solving generalized Stefan-Maxwell 
equations (ref. 8) involving nonzero m e and Maxwell’s relations. 

It is not obvious that the I, (1/p,) J, F„ term in equation (18) could indeed be neglected especially 
when J = 0, where electron and heavy particle diffusion velocities would be comparable. ( he 
I C, (u - u >F term in equation ( 1 5) is not negligible, since diffusion velocity of electrons for 
nonzero J can be large.) Nevertheless, evaluation of those terms does not represent additional 
complexity in multicomponent MHD, since friction coefficients a :j (or diffusion coefficients D (> ) and 
diffusion fluxes J, are obtained anyway, when species conservation equations are solved. 
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